Domain Decomposition Methods and High-Order Absorbing Boundary Conditions for the Numerical Simulation of the Time Dependent Schrödinger Equation with Ionization and Recombination by Intense Electric Field

This paper is devoted to the efficient computation of the Time Dependent Schrödinger Equation (TDSE) for quantum particles subject to intense electromagnetic fields including ionization and recombination of electrons with their parent ion. The proposed approach is based on a domain decomposition technique, allowing a fine computation of the wavefunction in the vicinity of the nuclei located in a domain Ω1 and a fast computation in a roughly meshed domain Ω2 far from the nuclei where the electrons are assumed free. The key ingredients in the method are i) well designed transmission boundary conditions on ∂Ω1 (resp. ∂Ω2) in order to estimate the part of the wavefunction “leaving” Domain Ω1 (resp. Ω2), ii) a Schwarz waveform relaxation algorithm to accurately reconstruct the solution. The developed method makes it possible for electrons to travel from one domain to another without loosing accuracy, when the frontier or the overlapping region between two domains is crossed by the wavefunction.